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O ■ The Wigner-Weisskopf-type model developed in [R. Alicki and F. Giraldi, J. Pliys. B 

cni: 

5_i . 44, 154020 (2011)] is applied to the biological process of energy transfer from a large 

O-i 
^ \ peripheral light harvesting antenna to the reaction center. This process is mediated 

•/^ ■ by the Fenna-Matthews-Olson (FMO) photosynthetic complex with a remarkably 

' (-h '! high efficiency. The proposed model provides a simple resonance mechanism of this 

phenomenon employing exciton coherent motion and described by analytical formu- 
las. A coupling to the vibrational environment is a necessary component of this 



e3 

qh| mechanism as well as a fine tuning of the FMO complex Hamiltonian. The role of 



the relatively strong coupling to the energy sink in achieving the resonance condition 
^ . and the absence of heating of the vibrational environment are emphasized. 
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I. INTRODUCTION 



The almost perfect energy transfer (APET) in photosynthesis phenomena, reaching the 
values over 95%, attracts attention of physicists in the recent yearsi. In particular, the ex- 
perimental observations of long-lasting quantum coherence in such systems^i^ suggest that 
the popular picture of the excitation moving stepwise from the highest to the lowest exci- 
tonic state, by means of a diffusive hopping, need not to be correct. A natural consequence 
of this observation is that the quantum theory of open systems should be applicable to these 
phenomena. Indeed, a number of papers were published where the interplay between quan- 
tum coherence and environmental decoherence was considered as a plausible explanation of 
APET. Those papers present a variety of theoretical approaches which differ by a choice 
of mathematical modeling of an environmental noise. In some papers Markovian models of 
Lindblad-Gorini-Kossakowski-Sudarshan are used^"™, in others non-Markovian effects has 
been also included using hierarchical equations of motioni^"— or certain "time-nonlocal" 
master equations^^. In all those models reaction center plays a passive role of a probability 
sink. Different regimes can be observed, depending on the choice of parameters, and various 
optimization schemes are studied. The overall picture shows a noise-assisted enhancement 
of energy transport accompanied by suppression of coherent oscillations, with the excep- 
tion of^^, where a small ensemble of molecular configurations yielding an efficient quantum 
transport (slightly reduced by Markovian dephasing), is found. 



In the present paper a completely different model is proposed - a resonant, coherent en- 
ergy exchange between the donor site and the acceptor one accompanied by an irreversible 
energy transfer to the sink (reaction center). The sink plays an active role modifying the 
energy of acceptor to achieve a resonance with respect to a donor. The applied mathematical 
formalism, based on the version of Wigner-Weisskopf model, developed in^^ is used. In con- 
trast to often very elaborate numerical analysis of master equations, or even the Schrodinger 
equation for the FMO complex and the bath, our approach employs simple analytical for- 
mulas which explicitly show a resonance APET mechanism and the necessity of fine-tuning 
of certain parameters. 



II. A MODEL OF "BOUNCING EXCITON" 

The structure of monomeric subunit of the FMO complex (see2^, Fig. 5) can be described 
as a cavity made of proteins containing the BChl pigments and suggests the following analogy 
proposed in^^. Consider an optical cavity which consists of two parabolic mirrors with a 
common symmetry axis and two identical 2-level atoms placed in the focuses of the mirrors 
with transition dipoles parallel to the axis. A simple physical intuition suggests that when 
one of the atoms is excited the emitted photon will be bouncing between atoms during its 
lifetime which depends on the cavity quality factor. If additionally the first atom (donor) is 
excited by a light source (say a laser) and the second (acceptor) is coupled to an energy sink 
we obtain a mechanism of energy transfer from a source to a sink mediated by the donor- 
photon-acceptor system. The similar mechanism can be responsible for the efficient energy 
transfer in photosynthesis phenomena. Two outer BChl pigments coupled to the antenna 
and the reaction center, respectively, are analogs of two atoms - donor and acceptor, and 
the other pigments together with the surrounding proteins form an environment which is 
the analog of electromagnetic field and supports wave-like traveling excitons - the analogs 
of photons. Under certain resonance conditions, which for the optical analog are guaranteed 
by identical frequencies of both atoms, the symmetry of the cavity and the symmetry of 
atoms positions and orientations, the "bouncing exciton" can transport energy from the 
antenna to the reaction center. The efficiency of the process depends on two parameters p 
and q. Here, p is the probability that the energy which reaches the acceptor in "one-shot" 
is dissipated to the sink and q is the probability of exciton recombination during the flight 
from the donor to the acceptor (or back). The simple analysis leads to the following formula 
for the efficiency of energy transfer after n bounces 

Using, the relation q ~ to/r where to — Ips is a typical time of flight and r ~ Ins is a typical 
exciton recombination time one obtains q ~ 10~^. In the next Sections we show that under 
reasonable conditions p ~ 0.5. Hence, using the above choice of parameters we obtain the 
following estimations for the "flve-steps" and asymptotic efficiencies 

ri{5) ~ 0.97, r/(oo) ~ 1 - g [l + ^ilz^] ~ 0.997 (2) 

which are sufficiently high to account for the observed in Nature values over 0.95. 



III. A MODEL OF EXCITONIC ENERGY TRANSPORT 

We consider a standard tight-binding model of energy transport in quantum networks 
in the single-exciton approximation which is vahd under the assumption that the exciton 
hfetime is much longer than any other relevant time scale in this model. The initial "skeleton 
system" consists of A^ + 1 sites with the corresponding single-exciton Hilbert space C^"^^ 
and the Hamiltonian matrix H = [hkf, k,l = 1,2, ...,N,N + 1]. The following assumptions, 
notation and terminology will be used: 

• The diagonal elements h^k are energies of the sites (denoted by Uk = h^k) while 
{hki; k ^ 1} are hopping amplitudes. 

• The sites " 1" , "2" , "N-l-1" correspond to a donor, an acceptor and a sink, respectively. 

• The direct coupling between the donor and the acceptor is usually small because of 
their spatial separation and similarly the sink is strongly coupled only to the acceptor, 
so we can put 

/ii2 = 0, /ij,^+i = for jV2 (3) 

what essentially simplifies the model. 

Under these conditions the Hamiltonian can be recasted into the form which reflects the 
open system character of the model 

H = u;i|l)(l| + uj2\2){2\ + Hb + {\l){gi\ + \2){g2\ + h.c.) . (4) 

Here a donor and an acceptor form a system with the Hamiltonian given by the first two 
terms of (jl]), the "skeleton bath" consists of the rest A^ — 1 sites including a sink with the 
Hamiltonian Hb- The last terms describe the system-bath coupling, |1) = [1,0, ...,0], |2) = 
[0,1,. ..,0], l^fi) = [0,0,hi3,...hiN,0], \g2) = [0,0,h23, ...h2N,h2,N+i] and Hb is a submatrix 
of H with indices k,l = 3,4, ...,N,N + 1, namely 
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Hb can be written in its spectral decomposition form 

7V+1 



Hb = 2_. ^o.\(y){oi\ 



(6) 



a=3 



where |A^ + 1) = [0, ...0, 1] and e^+i = ujy+i correspond to the sink whose energy should be 
the lowest one and well-separated from the others. 

However, this model with the low dimensional Hilbert space and hence the Hamiltonian 
with well-separated energy levels cannot describe irreversible energy transfer phenomena. 
We have to take into account the coupling to vibrational degrees of freedom. This interac- 
tion transforms the eigenstates |a) into resonances with finite spectral widths (decoherence 
rates) F^. The finite-dimensional Hilbert space spanned by {\a)} is replaced by an infinite- 
dimensional L^(M) and the Hamiltonian ([6]) by the continuous spectrum multiplication op- 
erator by u. We do not need to determine the exact form of the total Hamiltonian denoted 
by H' but the whole information will be included in the form of the relevant correlation 
functions^^. The system begins its evolution in the donor state |1) and after time t can be 
found in the acceptor state |2) with the probability 



Vuit) = \Au{t)\\ A2 = (2|e-^^*|1). 



(7) 



The APET takes place if for a certain time to the transfer probability Vi2{to) — 1- Then the 
energy must be transferred from the acceptor to the sink in an irreversible process which is 
fast in comparison to the lifetime of the acceptor state. 

The relevant correlation functions {gj\gj{t)) are now given by 

N+l 

G,,{t) = Y,M9,)\'e-'"''e-^-\ J = 1,2, t>0 (8) 



a=3 



where the presence of decoherence rates Fq, reflects the fact that the states \a) are not the 
eigenstates of the modified Hamiltonian H'^ but the unstable resonances embedded into 
practically continuous spectrum. 

The Laplace transforms of the correlation functions read 
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iu]t, 



e'^'Gjj {t) dt = jj{uj) + i 6j (u) (9) 



7,M = ^|(a|,,)r ; (10) 

Af+1 

^,M = ^|(a|,,)r _"^~^" (11) 

The final approximate formula for the optimal transfer probability based on the derivation 
presented in^^ which is valid under the condition that the bath-mediated interaction and the 
collective dissipation effects between the donor and the acceptor can be neglected, reads 



with the decay rates 



and the energy shifts 



Pi2(to) 
Here 
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(12) 



^^■(^) = ^] f^\-^ i ^2 ' ^■ = 1'2 (13) 

T^\uj-ujj-bj{uj)\ +7i(^)' 

are interpreted as normalized^^ moduli squares of the wave packets emitted to the bath from 
the donor (forwards in time) and from the acceptor (backwards in time), respectively. The 
oscillating term exp[26'(co'; to)] containing the relative phase of these two wave packets reduces 
the value of transition probability and will be discussed in the next Section. 

IV. THE ROLE OF SINK AND PHASE-FACTOR 

In order to give a more transparent picture of APET we approximate the functions fjiiS) 
in (fT3|) by Lorentzians 

TT (w - Wj) + 7^^ 

where 7j = 7j(wp and the renormalized frequencies satisfy the following self-consistent 
equation 

u:] = u,^bM^. (15) 



The APET condition means the perfect overlap of both Lorentzians what imphes 

uj\ = u\ = Wo, 7i = 72 = 7- (16) 

Notice that due to the condition /ii^at+i = (A^ + l|(7i) = the renormahzed frequency uj^ and 
the decay rate 71 of the donor do not depend on the sink parameters /i2,Ar+i = (A^ + 1|5'2), 
eAT+i and r^r+i. Therefore, the APET conditions fITB]) can be reahzed by fitting the sink 
parameters to modify properly the renormahzed frequency iij\ ^^"^ "^^ decay rate 72 of the 
acceptor. The resonance condition (fT6|) implies that the exciton energy at the donor is the 
same as at the acceptor and equal to /iwo what means the absence of heating of the vibrational 
environment. To achieve a rapid and irreversible energy transfer from the acceptor to the 
sink, the acceptor-sink coupling /i2,iv+i and the decoherence rate for the sink Fjv+i should be 
much larger than the decay rate 7. Moreover, the parameter PAf+i can be treated as a sum 
of two terms, Fat+i = P^^^^ + P^+i where P^+i corresponds to the decoherence phenomena 
while P7V+1 describes exciton decay in the sink. In fact one should add — iP^^^ to the sink 
energy e^z+i and then transform the initial Hamiltonian into a nonhermitian one to include 
this decay process. Such interpretation does not change the basic formula flT^ at resonance 
as the imaginary part of the phase Q{yi\ to), due to — iP^_,_^, appears around the sink energy 
cat+i which is assumed to be far away from the resonance value uq. On the other hand it 
provides an explanation of the energy transfer irreversibility. 

Finally, one has to discuss the influence of the oscillating factor which under the conditions 
(IT^ and fITB]) possesses a much simpler structure derived using the results o&- 

6{uj; to) = 2 arctan I I — {t{uj) — to) {uj — wo), (17) 

V 7 / 

where we use the parametrization of the phase-factor in terms of t{uj) which can be inter- 
preted as the propagation time of the exciton wave packet, with frequencies concentrated 
around u, between the donor and the acceptor. For a constant exciton velocity, i.e., t{uj) = r 
an optimal choice of to = r + l/7 gives the transition probability 4/e^ ~ 0.54. Simple numer- 
ical estimations show that for the choice t{uj) = t{u)o) + k,{u) — uq)'^ one can reach probability 
0.72. These relatively low probabilities are caused by the process of exciton back-scattering 
by the acceptor coupled to the sink. In order to explain the observed efficiencies higher 
than 0.95 we need the "bouncing exciton" model described above. In the following we shall 



concentrate on the values of the overlap integral 

T= j y/Mu)f2{u)du . (18) 

For example, if J-" ~ 0.7 and the presence of the phase flT7|) introduces another factor 0.7, 
one can estimate the single-shot transfer probability as p ~ 0.5 which is high enough to 
account for the observed phenomena. 

V. RESULTS FOR FMO COMPLEX 

We apply the presented theory to the case of the FMO complex of C. tepidum using the 
data of^"^. The model of FMO complex is given by the Hamiltonian of type ([5]) with N = 7 
containing 7x7 Hamiltonian submatrix taken from^. We adopt a standard choice of the 
BChl 1 and BChl 3 molecules as the donor and acceptor, respectively. Notice, that in the 
recent literaturei^'i^'2^'2^ a new model containing eight BChl pigments is applied, however, 
with a still large variance between Hamiltonian parameters found in the literature. As the 
aim of the present work is not to give the ultimate quantitative description, but rather to 
illustrate the new mechanism of energy transfer, we use the "old model". 

The free parameters of the model are randomly sampled within the following intervals: 
Us G [—500,0], Tj G [50,90], j = 1,...,8 (all values in cm~^). The limits put on the 
decay rates F^, j = 1,...,8, are consistent with the experimental widths of absorption 
spectral lines obtained at the temperature Tq = IIK^. The bound for wg is comparable 
to the largest energy variation in the FMO complex and the value of \h2%\ is practically 
unrestricted. Notice, that the value of ws must be negative, in order to "push up" acceptor 
energy to achieve resonance with a donor. 

Under these conditions the maximal value of the factor J-" ^ 0.75 and the corresponding 
resonance phenomenon at the point ujq ~ 150cm~^ is clearly visible (see Fig. [1^) with the 
width 7 ~ 30cm~^. The interesting feature is that the curves representing the shifted en- 
ergies for the donor and the acceptor are very close in a rather wide interval of frequencies 
[100,300] (Fig. Wp)- a similar behavior is observed for decay rates as well (Fig. [2^). It 
suggests that the Hamiltonian of FMO complex is well-tuned to make the resonance condi- 
tion less sensitive with respect to the parameters \h2s\ and Wg. The efficiency dependence 
on these parameters is illustrated by Fig. [3] and in Fig. [T]d. 

8 



The predictions concerning the sink parameters are the following: 

1. The ratio |ci;8|/|/i28| — 1-5. 

2. The value of |/i28| or — wg should reach maxiniuni available by the physical mechanisms 
in such systems. 

For our choice of the limit for cug the optimal value of acceptor-sink coupling |/i28| is about 
three times larger than the largest value of the dipole-dipole coupling between pigments 
in the FMO complex. This can be explained by the expected close contact between the 
acceptor and the reaction center—. 

To illustrate the temperature dependence of J-'(T) we assume that all decay rates are 
proportional to the environmental temperature^ and hence Tk{T) = TkiTo) x [T/To),k = 
1, . . . , 8. The results, with the optimal choice of rfc(To) as in Fig. [H are illustrated in Fig. 

SI 

VI. CONCLUDING REMARKS 

In contrast to the models based on Markovian or non-Markovian relaxation towards 
the lower energy levels in our model the exciton energy is not dissipated into vibrational 
environment of the FMO complex but is transfered in a resonant process to the acceptor 
and then transformed into other forms of energy in the reaction center. This implies the 
absence of heating of the vibrational degrees of freedom, which is considerably large in the 
standard models. Namely, in those models the difference between donor and acceptor energy 
is dissipated into environment before the exciton can reach the reaction center. This portion 
of energy ~ 200cm^^ ~ /c^ x 300i^ per exciton, can have biological relevance and, perhaps, 
could be used as a test for the proposed models. On the other hand, the vibrational degrees of 
freedom provide the necessary "quasi-continuous medium" for the energy propagation. The 
process of energy transfer is a coherent one and its efficiency is increased by the mechanism 
of "bouncing". Temperature dependence of the efficiency shows its rapid enhancement with 
increasing temperature and the relative stability in the vast region of high temperatures. 
The fine-tuning concerns mainly the Hamiltonian of the FMO complex what suggests an 
evolutionary natural selection mechanism. All these features seem to be consistent with the 
phenomenology of energy transfer in the FMO complex. 
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FIG. 1. The functions fi{uj) and y/f^JJJffoiiJj} for optimized parameters (in cm ^): r3= 59.6, r4= 
90.0, r5= 50.3, r6= 59.7, r7=89.7, r8=50.1, /i28=327 and ujs=-500 giving T ^ 0.75 (a) and with 
shifted /i28 = 160 resulting in T « 0.18 (b). 
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FIG. 2. The functions ^i{uj) (a) and 6i{u;) (b) for optimized parameters (see Fig. [T]) giving J-" ~ 0.75. 
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